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Abstract — In this paper we present an efficient Coupled-wave 
surface-impedance method for the analysis of extraordinary 
optical transmission (EOT) through single and stacked realistic 
metallic screens under normal and oblique incidence, including 
possible dielectric interlayers. The proposed theory is valid for 
the complete frequency range where EOT has been reported, 
, including microwaves and optics. Electromagnetic simulations 
validate the results of the model, which allows for a fast and 
accurate characterization of the analyzed structures. 

Index Terms — Extraordinary optical transmission (EOT), fish- 
' nets, surface impedance, mode-matching. 

I. Introduction 

THE phenomenon of Extraordinary Optical Transmission 
(EOT) through metallic screens with a periodic array 
of sub wavelength holes, first reported in [1], still is a topic 
of intensive research. The reader is referred to the excellent 
■ reviews by C. Genet et al. [2J, F. J. Garcia de Abajo O] 
and F. J. Garcia Vidal et al. [4J in order to have a complete 
overview of the topic. Coupling of surface plasmons through 
holes provided the first explanation for EOT iH, |I61. However, 
EOT has also been reported in cases where surface plasmons 
cannot exist, such as metals at microwave frequencies [7 1 and 
waveguide structures (81 . These phenomena show that EOT is 
a quite general effect that can be seen from many perspectives; 
plasmonics B|, @, diffraction models |10|, [HJ, 112J, [iJl, 
\ general surface- wave (or spoof-plasmon) analysis [jSJ, [fT4ll . 
and circuit analysis fTS], fT6l; an approach that has been 
promisingly extended to plasmonic structures [17|. Recently, 
some of the authors of this paper presented an analytical theory 
of EOT through perfect conducting screens, based on waveg- 
uide analysis |18|. This theory was later extended to screens 
made of realistic conductors and metals at optical frequencies 
1191 making use of the surface impedance concept, which 
is widely used in classical electromagnetism for analysing 
the skin effect in imperfect conductors [|2Ql . These theories 
provide a different and valid perspective on EOT but, due to 
the inherent simplicity of the models, they still have strong 
limitations: they are only valid for normal incidence and for 
holes smaller than - approximately - a quarter wavelength. 
These limitations arise from the approximations made in 
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the model, mainly from the "small hole approximation" first 
suggested by Gordon in [21]. 

EOT may find direct applications in photonic circuits [221 
and optical sensing ll23ll . EOT also is closely related to fishnet 
metamaterials El, EH, EH, EH. Complete fishnet analysis 
requires repeated computation of many stacked periodically 
perforated screens excited at oblique incidence [|27l . [|28l . 
Therefore, a fast and accurate method for the analysis of a 
single EOT screen is an unavoidable step for the complete 
characterization of fishnet metamaterials. In [|28l a mode- 
matching analysis of stacked perfect conducting EOT screens 
was shown to provide the necessary fast and accurate tool 
for this computation. However, direct application of this 
method to realistic metallic screens is cumbersome, due to the 
non-negligible penetration of fields inside metals at optical 
frequencies. In order to conveniently take into account this 
penetration without a substantial increase in complexity, the 
surface impedance concept offers an useful alternative [fT9ll . 

Coupled- wave method |29|, in general applicable also to 
non-uniform lines |30|, is considerably simplified in our case 
that only exhibits piecewise variation of cross-section. In 
such cases, the fields on the structure are simply expanded 
in the eigenmodes of the infinite waveguides matching the 
cross section at each point, except inside the metal, where 
we use approximate surface-impedance boundary conditions 
EH, lfT9ll in order to relate the electromagnetic fields on the 
interfaces. For the specific case of screens made of perfect 
conductors, the method reduces to standard mode-matching 
analysis f32| with n x E = boundary conditions on the 
metallic surfaces. 

Using scattering matrix formalism we will further extend 
the method to the case of stacked screens including dielectric 
substrates. This will allow us to study interesting phenomena 
such as anomalous EOT ll33l or "negative refraction" of 
Gaussian beams ll34ll . E3- 

II. Theory 

For simplicity, we will first consider a metallic screen with a 
periodic 2D array of square holes (the extension to rectangular 
holes and periodicities is straightforward). The geometry of 
the unit cell is shown in FigH] Oblique incidence of TE and 
TM waves, with tangential E ox H fields polarized along one 
of the main axes of the structure will be considered. Both 
cases are summarized in FiglTJb) and FigHJc). The transversal 
components of the electromagnetic fields at both sides of 
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Fig. 1. Front (a) and side views (b, c) of the unit cell of the structure with 
the three regions in which the fields are expanded. The incident waves in the 
cases of TE and TM polarization are shown in (b) and (c), respectively. Due 
to the polarization of the incident waves and due to periodicity, in (b) there 
are virtual electric walls in y = 0,±a/2 and in (c) there are virtual magnetic 
walls in a; = 0, ±a/2. The rest of lateral boundaries are periodic boundaries. 



the screen (regions 1 and 3) and inside the hole (region 2) 
are expanded in terms of Bloch and waveguide modes, and 
evaluated at the input and output surfaces of the screen (see 
Eqs. (ll4H17l) in the Appendix). The next step is to determine the 
coefficients in (I14H17I) . For this purpose we impose a surface 
impedance condition on the metallic interfaces at z = ±t/2 
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In ([3]), k^J^lf is the longitudinal component (i.e. parallel to 
the z axis) of the propagation constant of the "refracted wave" 
inside the metal, and Y 
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is the wave admittance of this 
"refracted wave". In this context, the "refracted wave" is 
defined as the wave that would be refracted at the air-metal 
interface if there were no holes in the screen, and the wave 
admittance is defined as the ratio between the transverse (in 
the X — y plane) components of the H and E fields of this 
refracted wave. For an incident TE wave as depicted in Fig. 
Kb) 
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where ^||^^' ^11^^ ^11^^ vectors containing 

the transversal components of the fields in regions 1 and 3, 
and Z is the surface impedance matrix corresponding to an 
infinite slab (without holes) for a given incident plane wave 
tin . When applying ([T]) we are neglecting the perturbation 
associated to the presence of the holes. Another assumption is 
that all the refracted modes can be connected by approximately 
the same surface impedance, so that ([B is satisfied by the total 
fields. This last approximation is valid as long as the transverse 
wavenumbers of the different modes (/ci^n, ky]}n) in ([T4l) and 
([T6b ) are much smaller than the TEM wavenumber inside the 
metal. This approximation will fail near the plasma frequency, 
which happens only for very small periodicities or for very big 
angles of incidence (see Section III. A). 

Matrix Z is symmetric, and can be diagonalized, so that 
the following combinations of the tangential components of 
the fields vanish at the screen interfaces [ 191 
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and for an incident TM wave as depicted in Fig. [He) 
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where k = w^emetMo and = i/emet/Mo are the 

propagation constant and the intrinsic admittance of the metal 
with relative permittivity Emet- 

In the area of the holes, continuity of transverse electro- 
magnetic fields is imposed, 
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for|x|,|y| <6/2. 



Equations Q and (O provide a complete set of boundary 
conditions. However, they are not appropriate for the direct 
computation of the coefficients in (I14H17I >. It is better to define 
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and combine ^ with (|6l|7]) to obtain 
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forO < \xl \y\ < a/2. 

These equations provide a complete set of four boundary 
conditions in tlie range b/2 < \x\,\y\ < a/2. We still need 
four more linearly independent equations in order to have a 
complete set of boundary conditions in the region |y| < 
b/2. For this purpose we choose 

Hl'\z = -t/2) = Hi;'\z = -t/2) 

hI^^ {z = t/2) = hI^^ {z = t/2) (9) 

for|x|,|^| <V2. 

Equations dg]) and (O are a complete set of boundary 
conditions, which provide a sparse system of equations for the 
mode coefficients (see Appendix). The specific case of perfect 
conducting screens is obtained by simply taking Z^i = Z^2 = 
instead of © . 

III. Results 

A. Single EOT screens 

In Figs. [2ll7] the transmission coefficients obtained with the 
reported model are compared with electromagnetic simulations 
using CST Microwave Studio. In order to obtain meaningful re- 
sults, the resolution of the highest mode inside the holes (P, Q) 
must be similar to that of the input and output regions (N^ M), 
i.e. [N,M] ^ [P,Q]{a/b). We employed P = Q = 2, which 
was enough to obtain accurate results compared with full wave 



0.8 

^0.6 



go.4 



g 0.2 
H 



= Q\ 




\a/4 






^ 


::Jj 


. -|[... 
iy-i-J 




^ m 











0.2 



0.1 









a/8 


///' \ \ 









0.995 1 

Ww.Ol 

(a) 



0.9998 1 

ffw, 01 

(c) 



Fig. 2. Transmission through an array of square holes in a copper screen 
(cr = 59.6 X 10^ S/m) at normal incidence and different sizes of the holes. 
Periodicity is a = 300/im and thickness of the screen is t = a/20. Wood's 
anomaly frequency is /w,oi ~ 1 THz. Continuous lines correspond to the 
mode matching model, dashed lines to the CST simulations and dotted-dashed 
lines to our previous model |19|. In (a) a general picture is shown; in (b) and 
(c) the details of the peaks for b = a/2 and b = a/8 are shown. 
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Fig. 3. Transmission through an array of square holes in a silver screen 
(ojp = 27r X 2175 THz and fc = 4.35 THz) at normal incidence and different 
sizes of the holes. Periodicity is a = 1/im and thickness of the screen is t = 
a/ 20. Wood's anomaly frequency is /w,oi ~ 299.79 THz. Continuous lines 
correspond to the mode matching model, dashed lines to the CST simulations 
and dotted-dashed lines to our previous model |T9]. 



electromagnetic simulations. CPU time per frequency point 
was ~ 0.5 s vs ~ 4 min with the electromagnetic solver. 
If more modes are used, computational time grows very fast 
(with P^), making the theory unpractical. In addition, there 
appear numerical errors due to very small amplitudes of high 
evanescent modes. 

In Figs. [21 131 and [5] the metallic screen is modeled by a finite 
conductivity a = 59.6 x 10^ S/m (corresponding to copper), 
and the electric permittivity is given by 

In Figs. O |6l and [71 the metallic screen is modeled by the 
Drude-Lorentz permittivity 
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with the plasma frequency cjp = 27r x 2175 THz corresponding 
to silver I136L The frequency of collision corresponding to 
silver was corrected by the factor (1 + /eff)/(2t) (where /eff 
is the mean free path in bulk silver) in order to take into 
account the finite thickness of the screen |19|. That is, we 
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Fig. 4. Transmission through an array of square holes in a copper screen 
(cr = 59.6 X 10^ S/m) at oblique incidence (TE incident wave). Periodicity 
is a = 300/im, hole size is 6 = a/4, and thickness of the screen is t = 
a/ 20. Continuous lines correspond to the mode matching model and dashed 
lines to the CST simulations. In (a), the peaks correspond to the divergence 
of the scattered TMo,i mode and the Wood's frequencies range from 1.00 
THz for 6* = 5° to 7.46 THz for = 60°. In (b), the peaks correspond 
to the divergence of the scattered TM_i,i mode and the Wood's anomaly 
frequencies range from 1.33 THz for 6 = 5° to 1.00 THz for 6 = 80°. 
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Fig. 5. Transmission through an array of square holes in a copper screen 
(cr = 59.6 X 10^ S/m) at oblique incidence (TM incident wave). Periodicity 
is a = 300/im, hole size is 6 = a/4 and thickness of the screen is t = a/20. 
Continuous lines correspond to the mode matching model and dashed lines 
to the CST simulations. In (a), the peaks correspond to the divergence of the 
scattered TMo,-i mode and the Wood's frequencies range from 0.92 THz 
for ^ = 5° to' 0.50 THz for = 80°. In (b), the peaks correspond to the 
divergence of the scattered TMo,i mode and the Wood's anomaly frequencies 
range from 1.09 THz for 6> = 5° to 1.35 THz for = 15°. In (c), the peaks 
correspond to the divergence of the scattered TMi,_i mode and the Wood's 
anomaly frequencies range from 1.27 THz for 6 = 10° to 1.17 THz for 
e = 20°. 



have chosen /c = (1 + ^eff)/(2t) x /c,Ag ^ 1.26 x 4.35 
THz. Moreover, in order to take into account the effect of 
field penetration through the walls of the holes, the size of the 
holes was increased to an effective width 6eff = b-\-26, where 
6 is the penetration depth of the fields inside the metal. For 
lossy metals at microwave or terahertz frequencies this field 
penetration is negligible (S « b), however for solid plasmas 
it has a relevant effect. The penetration depth at optical 
frequencies is evaluated sls S = Ap/(27r), where Ap is the 
plasma wavelength of the metal. This correction has effects on 
the frequency of the peaks, but has no effect on the estimation 
of losses. We will assume that losses inside the holes can be 
neglected with regard to losses in the external metal interfaces. 
The accuracy of this assumption is checked in the following, 
by comparison with electromagnetic simulations. 

In the Figures, the frequencies are normalized to those 



corresponding to Wood's anomalies (/w,nm) for each angle 
of incidence. These Wood's frequencies are defined as those 
frequencies at which the admittances of the different TMn,m 
modes scattered by the screen diverge. Therefore, they corre- 
spond to the solutions of the equation 

flnmiO) = + U,nmSm{0)y + (^)' (12) 

for an incident TE wave, and 

flnmiO)=(^Y+(^ + U,nmMO)y (13) 

for an incident TM wave. 

Figs. [2l|3] correspond to the case of normal incidence and 
different geometries. The coupled- wave model presented here 
is compared with the "small holes" model reported in [ 19], and 
with the results obtained using a conventional electromagnetic 
solver, which were used as reference. As expected, the "small 
hole" model fails for screens with holes larger than 6 = a/4 
(see Fig. O. The proposed new model provides good results 
even for holes as large as 6 = a/ 2. For larger holes it starts to 
be inaccurate. However, in those cases the transmission peak is 
due to the resonance of the hole and can not be properly called 
"extraordinary". As was mentioned earlier, these results were 
obtained with computation times several order of magnitude 
smaller than those obtained using conventional electromag- 
netic solvers. 

In the remaining figures - Figs. |4] to |7] - oblique incidence 
of TE and TM waves over copper and silver screens operating 
at THz and optical frequencies, respectively, is considered. A 
quite good agreement is found between the reported model 
and the simulations for any angle of incidence, periodicity 
and hole size. Only in cases of very high angles of incidence 
and frequencies high enough to provide relatively small values 
of permittivity (Ismetl ~ 1 or 10), the results of the model 
deviate significantly from those of the simulations. These 
cases, however, are of small practical interest, because the 
background transmission is high and the height of the peaks 
is small. 

Since the Wood's frequency values (ll2H13t are different 
for the different angles of incidence, and the frequencies 
are normalized to Wood's values in the figures, the actual 
frequencies corresponding to the peaks shown in these figures 
cover a quite wide range of values (see captions). For an 
incident TM wave impinging at an angle d = 11.54^ on the 
screen, the Wood's frequency corresponding to the divergence 
of the scattered TMo,i and TMi^_i modes is the same. For 
higher angles of incidence, the peaks corresponding to the 
divergence of the TMo,i mode become weaker (Figs. Ob) and 
Ob)), and those corresponding to the divergence of the TMi^_i 
become stronger (Figs. EJc) andQc)). A similar effect occurs 
for incident TE waves at i9 = 26.57^ and TMo,i and TM_i,i 
modes. 

The surface impedance condition ([T]) used along this analy- 
sis is also properly taking into account the tunneling through 
the metal, when the screen thickness becomes of the same 
order or smaller than the penetration depth. This effect is 
studied in Fig. [8] where EOT through silver screens of dif- 
ferent thicknesses is studied. As expected, tunneling through 



5 



0.6 



0.4 



a 

g 0.2 
H 





















lU 

i \ 








0.95 

.//?;v.o. 

(a) 



Fig. 6. Transmission through an array of square holes in a silver screen 
(cjp = 27r X 2175 THz and fc = 4.35 THz) at oblique incidence (TE 
incident wave). Periodicity is a = l/xm, hole size is b = a/ 4: and thickness 
of the screen is f = a/20. Continuous lines correspond to the mode matching 
model and dashed lines to the CST simulations. In (a), the peaks correspond 
to the divergence of the scattered TMo,i mode and the Wood's frequencies 
range from 300.93 THz for 6* = 5° to 318.98 THz for 6 = 20°. In (b), the 
peaks correspond to the divergence of the scattered TM_i,i mode and the 
Wood's anomaly frequencies range from 400.07 THz for ^ = 5° to 315.08 
THz for 6* = 40°. 
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Fig. 7. Transmission through an array of square holes in a silver screen 
(ujp = 27T X 2175 THz and fc = 4.35 THz) at oblique incidence (TM 
incident wave). Periodicity is a = Ifim, hole size is b = a/4 and thickness 
of the screen is t = a/20. Continuous lines correspond to the mode matching 
model and dashed lines to the CST simulations. In (a), the peaks correspond 
to the divergence of the scattered TMo,-i mode and the Wood's frequencies 
range from 275.51 THz for 6> = 5° to 151.10 THz for = 80°. In (b), 
the peaks correspond to the divergence of the scattered TMo,i mode and 
the Wood's anomaly frequencies range from 328.39 THz for ^ = 5° to 
455.68 THz for = 20°. In (c), the peaks correspond to the divergence of 
the scattered TMi^_i mode and the Wood's anomaly frequencies range from 
400.07 THz for = 5° to 299.79 THz for = 80°. 



the metal results in increased background transmission and 
smoother transmission peaks. The agreement of our results 
with the electromagnetic simulations is very good, with no 
loss of accuracy for thinner screens. 

B. Dielectric loaded screens 

Realistic EOT screens must be supported in most cases 
by a dielectric board. The presence of the dielectrics also 
introduces new effects such as anomalous EOT B3I . The 
analysis presented in Section HIl can be easily extended to these 
cases using the standard scattering matrix method. This implies 
to consider the presence of a complete set of electromagnetic 
modes impinging from z ^ ±oo in Fig. [T] and to compute the 
complete scattering matrix between the interfaces at z = ±t/2 
for all these modes. Once this matrix is computed, we only 
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Fig. 8. Transmission through an array of square holes in a silver screen 
{ujp = 27T X 2175 THz and fc = 4.35 THz) at normal incidence and 
different thicknesses of the screen. Periodicity is a = 1/im and hole size 
is b = a/4. Wood's anomaly frequency is /w,oi ~ 299.79 THz. Continuous 
lines correspond to the mode matching model and dashed lines to the CST 
simulations. 




Fig. 9. Transmission through an array of square holes in an aluminium layer 
(ujp = 27T X 3570 THz and fc = 54.11 THz) deposited on a dielectric slab 
(£d = 2.25 and tan 5 = 0.001). The incident plane wave is a TEM wave 
with the electric field polarized either along the short (the so called anomalous 
extraordinary transmission [33 ) and long (the so called regular extraordinary 
transmission |33 |) periodicities. Periodicities are 47.5 x 113/im, hole size 
is 35.2/im, thickness of the metallic layer is 0.5/im and thickness of the 
dielectric slab is 20/im. Continuous lines correspond to the mode matching 
model and dashed lines to the CST simulations. 



need to combine this matrix with the propagation matrices 
through the dielectric board(s), and the scattering matrices of 
the air-dielectric interface(s). This procedure is cumbersome 
but straightforward, therefore, it will not be shown explicitly. 
At this point it is worth noting that the cascading of scattering 
matrixes is not as simple as the cascading of transmission 
matrixes [32], but gives important advantage of the numerical 
stability with respect to ill-conditioned transmission matrixes 
which use both, the negative and the positive exponentials. 

In Fig. [9] EOT through a dielectric loaded metallic screen 
similar to that studied in ll33l (Fig. 3. a) is analyzed. In 
particular, we used the same periodicity and substitute circular 
holes by square holes of the same cross section, obtaining 
almost identical results, with the anomalous and "regular" EOT 
peaks clearly present. 

C. Refraction of Gaussian beams 

As it is well known, laws of refraction in material slabs 
can be determined by analyzing the behavior of light beams 
incident at an oblique angle over the slab. For very large 
(positive or negative) index of refraction, the beam is refracted 
along the normal to the slab interface, and the output beam 
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Fig. 10. Transmission through a single (a) and four stacked (b) copper screens 
(cr = 59.6 X 10^ S/m) at oblique incidence (TM incident wave). In (a) the 
metallic screen is sandwiched between two dielectric slabs (e^ = 2 A3) and 
in (b) the four metallic screens are placed among a total of five dielectric 
slabs. Periodicities are 1.5 x 3.4 mm, size of the square holes is 1.1 mm, 
thickness of the metallic layer is 35/im and thickness of each of the dielectric 
slabs is 0.49 mm. Continuous lines correspond to the mode matching model 
and dashed lines to the CST simulations. 

emerges from the opposite side of the slab with some lateral 
shift which depends on the slab thickness. Assuming this shift 
as the reference, positive or negative shifts with regard to this 
reference are the signature of positive or negative index of 
refraction. In [[34l and [|27l the behavior of electromagnetic 
Gaussian beams impinging at an oblique angle on single and 
stacked EOT screens was studied. It was found that in all 
cases TM beams experience a "negative" shift, whereas TE 
beams always experience positive or no shift, which was 
consistent with the presence of negative refraction in fishnet 
metamaterials |34|. 

The scattering matrix scheme in Section IIII-BI can be also 
applied to stacks of EOT screens. Once the transmission 
coefficients for incident plane waves are known, it is a rela- 
tively straightforward task to analyze the behavior of incident 
Gaussian beams, which can be decomposed into plane waves. 
In Fig.[TOlwe show the transmission coefficient for TM plane 
waves impinging at several angles over dielectric loaded single 
and stacked EOT screens, with an structure very similar to 
the structures studied in [|34l . As can be seen in the figure, 
our results show a very good agreement with the electro- 
magnetic simulations, with much smaller computation times. 
Once the position of the peaks have been determined, we are 
ready for the computation of the lateral shift of an incident 
Gaussian beam at the frequency of these transmission peaks. 
The profiles of the input and output TM beams impinging 
at different angles over the single and stacked structures are 
shown in Fig. [TT] Each output beam profile is centered on the 
aforementioned reference point, corresponding to an infinite 
"effective refractive index" for the structure. As can be seen, 
a negative shift is observed for all output beams, in agreement 
with the experimental results shown in [34 | and [27|. Similar 
computations have been carried out for TE beams (not shown), 
confirming only small positive shift. 

IV. Conclusion 

An analytical method has been provided for the analysis 
of extraordinary optical transmission (EOT) through realistic 
metallic screens with a periodic array of holes. Our method is 




Fig. 11. Electric field profiles of a ID Gaussian beam impinging at different 
angles of incidence and passing through the same structures as in FigfTol 
For each angle of incidence, the fields are calculated at the frequency of the 
maximum transmission. For the single structure (a) these frequencies are 61.2, 
60.7, 59.6, 56.3 and 52.7 THz for angles of incidence 0, 5, 10, 20 and 30° 
respectively. For the four-stacked-screens structure (b) these frequencies are 
52.1, 51.1, 49.2, 45.3 and 42.1 THz for angles of incidence 0, 5, 10, 20 and 
30° respectively. Deviation of the beam to the left is higher for increasing 
angles of incidence. The left inset shows the polarization of the beam relative 
to the screen(s) as well as the planes where the field profiles are depicted. 
The input and output planes are placed 100 mm away from the origin of the 
cartesian axes. 



based on the coupled-wave analysis with approximate surface 
impedance boundary conditions at the metallic interfaces. This 
method allows for a fast and accurate characterization of EOT 
for any angles of incidence, hole sizes, thicknesses and all 
frequencies where EOT has been reported (from microwaves 
to optics). Due to its computational efficiency, the method can 
be applied to the analysis of dielectric loaded and stacked 
EOT screens, including complicated sources, such as Gaus- 
sian beams. The method has been tested by comparing the 
results with those coming from commercial general purpose 
electromagnetic solvers (with computation times several orders 
of magnitude smaller) and experimental work on anomalous 
EOT and negative shift of Gaussian beams through single 
and stacked EOT screens. A very good agreement has been 
obtained in all cases. Future work will be aimed to the com- 
plete characterization of these and other exciting phenomena 
in single and stacked EOT screens and fishnet metamaterials. 



Appendix A 



The transversal components of the electromagnetic fields at 
both sides of the screen (regions 1 and 3) and inside the hole 
(region 2) can be expanded in terms of Bloch and waveguide 
modes and evaluated in the input and output surfaces of the 
screen. 

In case of an incident TE wave with the electric field along 
the y direction, the magnetic field along x and z and the 
propagation vector along x and z, the transversal electric field 
components in the input and output surfaces of the screen are 
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where i?™, i?™, T™ and T™ are the ampHtudes of the 
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different modes scattered in the screen, k^n = kx o -\ 

2m7r ' ' ^ 

and ky]}n = are the transversal components of the wave 

vectors of the different modes, and kx^o = A:osin(^) with 
being the angle of incidence. The transversal electric fields 
inside the holes in this case are 
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where S*™^ and 5*™^ are the amplitudes of the different 
modes inside the holes propagating either to the left (~) or 

to the right (+), and kl^l = ki^l = and ki^lq = 
' ' 

^h^o ~ {^x'}p^ — {k^}q^ are the components of the wave 
vectors of the different modes, e^^ is the relative permittivity 
of the material filling the holes. In all the simulations in this 



(2)/ 



U(2). 



paper = 1. Note that for this particular polarization there 
are virtual electric walls in^ = 0,=ba/2. 

In case of an incident TM wave with the electric field along 
the y and z direction, the magnetic field along x and the 
propagation vector along y and z, the transversal electric field 
components in the input and output surfaces of the screen are 
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transversal components of the wave vectors of tfie different 
modes scattered in the screen and ky^o = kosm{0), with 
being the angle of incidence. The transversal electric fields 
inside the holes in this case are 
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'X,p 



' — ~r~ and kz^pq — 



)J ShkQ — ^4^p^ — ^4^g^ are the components of the wave 
vectors of the different modes inside the holes. Note that for 
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this particular polarization there are virtual magnetic walls in 

X 0,±a/2. 

Expressions for magnetic fields Hx are obtained by just 
multiplying each mode in Eqs. (11411171) by its corresponding 
wave admittance. Expansions in Eqs. ([T5]) and ([TTl) correspond 
to a square waveguide with PEC walls, so losses inside the 
holes are neglected. 

Boundary conditions ([5]) and (O can be multiplied by a 
convenient set of orthogonal functions and transformed in 
integral boundary conditions over the screen and in the region 
of the holes. In the cross section of the unit cell, the following 
integral boundary conditions are satisfied 



conditions 



7^(3) , 7^(1) y 



JrLx — -tl: 



(1)11 . 



E^^ + — Zsi — H^^'^ I fnmdS — 



// { ^- 

J J wg 

-IIS ^ 



(3) 7^(1) y 



(2) T^{2) ry 





- J^x 







-m 



77'(3) I 771(1) I r7 



E^^ + E^^ -\- Zsi 



rrw; Try 

Uy - Hy 



r(3) 

y 



Uy - Uy 



I fnmdS 
} fnmdS = 

I QnmdS 
} 9nmdS = 



jj H^^ ■ Upq{x,y)dS = J J H^^ ■ Upq{x,y)dS 



Hi, ■ Wpq{x,y)dS = Hy ■ Wpq{x,y)dS 

h J Jh 
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(21) 



Hy^'^ • Wpq{x,y)dS ■■ 



J J H^^^ • Wpq{x,y)dS 



where Upq and Wpq are orthogonal functions in the area of the 
unit cell. For incident TE waves with the electric field along 

y 

Upq{x, y) = sm{kx,pix + 6/2)) cos{ky^qy), 

Wpq{x, y) = cos{kx,p{x + 6/2)) sm{ky^qy), (22) 

p=l,2,..,P; ^ = 1,2,..,Q; 
and for incident TM waves with the magnetic field along x 

Upq{x, y) = cos{kx,px) cos{ky^q{y + 6/2)), 



(x, y) = cos{kx,px) cos{ky^q{y + 6/2)), (23) 



p = l,2,..,P; g = l,2,..,Q. 

After substitution of (11411171) into ([T8]) and (|2T]), these inte- 
grals are solved analytically, providing an sparse system of 
equations for the coefficients in (I14H17I) . 



J J wc 



^s2 



Jwg " 

-IIS 



rjx — -t^x -r ^s2 



H^^ + I QnmdS = 0, 

(18) 

where subindex wg and h stand for the waveguide and 
hole sections respectively, and fnm and gnm are orthogonal 
functions in the area of the unit cell. For incident TE waves 
with the electric field along y 

fnm{x, y) = ex.p{-ik^]lx) cos(4!m^), 
gnm{x, y) = exp{-iki]lx) sm{k^]lny), (19) 
n = -A/',...,0,...,A/'; m = l,2,...,M; 
and for incident TM waves with the magnetic field along x 
fnm{x,y) = cos{k^]lx)ex.p{-ik^]}ny), 
gnm{x, y) = sm{k^]lx) ex.p{-ik^]}ny), (20) 
n = l,2,..,Ar; m = -M,...,0,...,M. 

In the area of the holes, continuity of transverse electromag- 
netic fields is imposed with the following integral boundary 
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